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Lecture Notes on Fluid Dynamics 

(1.63J/2.21J) 
by Chiang C. Mei, 2002 

CHAPTER 6. 

SEEPAGE AND THERMAL EFFECTS 
IN POROUS MEDIA 

6-ldarcy-EM.tex 

Applications : Groundwater flow and transport, building insulation, energy storage and re- 
covery, geothermal reservoirs, nuclear waste disposal,etc. 



6.1 Empirical basis of Darcy's law for seepage flow 

[References] : 

Polubarinova-Kochina: Theory of Groundwater Movement, Princeton University Press 
The basis of Darcy's law for a non-deformable medium is the one dimensional experiment 
by Darcy, see Figure 6.1.1. 




Figure 6.1.1: Darcy's experiment for seepage flow 
The discharge through the tube is measured to be 



Q = KA*^1 (6.1.1) 



2 




Figure 6.1.2: A one dimensional model of porous medium 



where 

P P 
Hi = h yi H 2 = h 2/2 

pg pg 

and K is an empircal coefficient called the hydraulic conductivity. The effective (seepage, or 
filtration) velocity is defined to be the discharge per unit gross area of the porous medium 

Q T ,H 2 - Hi 



Let us define the potential to be 



= -k \y g +y) = ~ KH (6- 1 - 3 ) 



<j>i = —K [ — + Vi) = -KHi, ^ = 1, 2, (6.1.4) 



then 

4>i = -K ( 

.pg 

In the limit of As — > 0, we get 

u = - k ^ = js (6 - L5) 

This is Darcy's empirical law relating the seepage velocity to the hydraulic heads, both are 
macro-scale averaged quantities. 

What affect the conductivity? Let the porosity n be defined as the percentage of pore 
volume in the gross volume V . If V s is the volume occupied by solid grains in V, then 

V - V. 

n=~ y± (6.1.6) 

If the pores are satuarated with fluid then Vf — V — V a and n = Vf/V. 

Consider the idealized porous medium consisting of parallel tubes, Figure 6.1.2. In a 
cross-section of area A, the net area of pores is nA. The net averaged velocity is 

Q U (p. i n\ 

Up = -A = n ( °- 7) 
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For a laminar flow through a circular tube, the discharge is 

(Homework). The averaged velocity in the pore (tube) is 

Q R 2 g (&P/P9 + Ay) 
Up irR 2 8is As 

The seepage velocity through the matrix is 



nR 2 g (Ap/pg + Ay) 

nup = 7 6.1.8 

8u As y J 



Therefore, the hydraulic conductivity is 



K = ^L (6.1.9) 

with the dimension 

[K] = ^ (6.1.10) 

The real pores are, of course, geometrically more complex, but the preceding formula indi- 
cates that K is small for small pores and for high viscosity. 
Darcy's law is also often expressed in the form, 

k 

u = — (p + pgy) (6.1.11) 

where k is called the (intrinsic) permeability. Clearly 

k = ^- (6.1.12) 
99 

For the tubular model we have 

k = ^f (6.1.13) 

which is independent of viscosity. 

From experiments, Kozeny and Carman proposed the following (empirical) formula 

k = cd 2 - (6.1.14) 

(1 — n) z 

where c = 0.1 ~ 0.8, and d is the effective pore diameter defined as the ratio of the volume 
of solids to the wetted area in the gross volume. 
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Materials 


Hydraulic 


Conductivity 


Clays 


< 10 _ 


9 




Sandy clays 


10- 9 - 


- 10' 


-8 


Peat 


10-9 - 


- 10" 


-7 


Silt 


io- 8 - 


- 10" 


-7 


Very fine sands 


io- 6 - 


- 10" 


-5 


Fine sands 


io- 5 - 


- 10" 


-4 


Coarse sands 


io- 4 - 


- 10" 


-3 


Sand with gravel 


10- 3 - 


- 10" 


-2 


Gravels 


> io- 
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fT(m/sec) 



Table 6.1: The order of magnitude of the conductivity of natural soils 



Three-dimensional Darcy law; As a generalization 

S , = = J^£±M (6.1.15) 

OXj \l OXj 

where denotes the conductivity tensor and kij the permeability tensor. For an macro- 
scopically isotropic material 

K i:j = KSij, = k5ij (6.1.16) 



Continuity requires that 
Thus 

in general, and 



V-m = (6.1.17) 



" r K,^)= (6.1.18) 



dxj V dx. 



V-KV(f) = (6.1.19) 
for isotropic media. If further, the material is homogeneous : K = constant, then 

V 2 = O (6.1.20) 

Hence, for a nondeformable isotropic and homogeneous porous medium, the flow is potential. 

In most soils the pore flow is usually laminar. Take the typical values : u = 0.25cm/ sec, d ^ 
0.4 mm , then 

uD 0.01cm 2 1 sec 

Re= = f. = 0.1. 

v 0.01 cm 1 j sec 

It is known empirically that for Re < 1 ~ 15, the flow is usually laminar. 
Boundary conditions of a typical seepage problem. Consder an earth dam: 



5 



On the soil water interface y = Hi(x) AB: 

P = Pa + P9 (Hi - y) 

(f) = -K ( — + y ) = -K ( — + E x - y + y ) = -K ( — + H 1 ) = constant 
\P9 J \P9 J \P9 ) 

Therefore, 

4> = constant, y = Hi(x). (6.1.21) 

On the phreatic surface AE, where y = Y(x) is unknown a priori. The dynamic condition is 
p = constant. Therefore, 

+ KY(x) = constant, y = Y(x) (6.1.22) 
In addition, we have the kinematic condition, that y = Y{x) is a streamline : 

ip = constant, y = Y(x) (6.1.23) 

On the seepage surface ED: 

p = <fi + Ky = constant = p a , y = H 2 (x) (6.1.24) 

where H 2 (x) is known. Note that just inside the soil, the fluid velocity is not tangential to 
ED. 

On the impervious boundary (rock or saturated fine clay, etc.), BF: 

ip = constant (6.1.25) 

Because the phreatic surface is unknown, the boundary-value problem in the earth dam 
is highly nonlinear and difficult. 

The prediction of the phreatic surface is useful for the design of the dam thickness and 
of the size and location of the drainage ditch. To find the stresses in the dam and on the 
foundation, one should consider the deformation of soil; this requires soil mechanics. 
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Figure 6.1.3: A one dimensional model of porous medium 



